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Abstract 

Within the RSJ model we performed an analytical and numerical investigation 
of SQUID cells consisting of two Josephson junctions shunted by an extremely small 
inductance leading to strong coupling of the elements. Contrary to the well-known 
behavior of cells shunted by a high inductance voltage phases of the junctions are 
locked with very small phase difference for almost all values of external flux. Only 
for external flux in the vicinity of half a flux quantum phase difference rises rapidly 

tO 7T. 
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There has been growing interest in the development of complex Josephson junction 
arrays during last years being caused by possible applications in the field of satellite com- 
munication systems, astronomical observations or construction of supercomputer chips 
[1-4]. Despite some promising experimental results with 2D arrays [5-9] detailed theoret- 
ical studies are still missing. In addition, nearly all theoretical investigations on coupled 
Josephson junctions in general are based on weakly coupled elements. Otherwise, present 
day and near future technology will allow construction of arrays with elements strongly 
coupled by inductances around 3pH and smaller. 

Thus before attacking more complex problems it is useful to perform theoretical stud- 
ies of very simple arrays of strongly coupled Josephson junctions. One of the simplest 
arrays at all is a completely superconducting SQUID cell consisting of two identical 
Josephson junctions with parallel current support inductively shunted by a very small 
(L <C 0o/2vrJc; O : flux quantum, Iq : critical current) inductance. Phase-locking 
in this model has been studied thoroughly in the weak-coupling limit several years ago 
]1|, [3], [Hj, but any detailed theoretical investigations on strong inductive coupling seem 
to be missing so far. The results of our work show that the behavior of strongly cou- 
pled Josephson junctions, especially under the influence of external magnetic fields 0, is 
qualitatively different from what has been known for weakly coupled elements. 

The circuit to be considered here is shown in Fig. 1. 

Both junctions, to be considered as being identical and having negligible capacitance, 
are described within the RSJ model. The behavior for weak coupling (L ^> 0o/27ric) 
can be summarized as follows: (i) In any case there exists a stable phase-locking regime, 
(ii) Without an external magnetic field the difference between Josephson phases on both 
junctions vanishes in the mean. In this case there is no stationary current flow through 
the shunt, (iii) Any external flux leads to a non- vanishing difference of Josephson phases 
of the elements causing a phase difference 5 between voltages. The mean value of this 
voltage phase difference 5 can be determined by the so-called reduced equation [1,10] 




i Q = I /I c , V9 = 2vr0/0 O - (2) 

Thus, for usual operation regimes with i « 1.5 the dependence is nearly linear (cf. 
Fig. 2). 

This may cause problems in larger arrays; if parasitic magnetic fluxes cannot be sup- 
pressed radiation output will be small because of cancellations between contributions of 
different cells. 

Here we consider the behavior for strong inductive coupling, L <C 0o/27rlc. For 
vanishing McCumber parameter j3c the circuit can be described by the RSJ equations 

4>i + sin0i = i - r 1 ^! - 02 + p), (3) 
02 + sin 02 = i + I' 1 (01 - 02 + V 9 ) (4) 
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where 

I = 2nI c L/(j) « 1 (5) 

is the dimensionless inductive coupling and the dot denotes the derivative with respect to 
the normalized time 

s = ^R N I c t. (6) 

For obtaining approximate solutions in the strong coupling limit one introduces new 
variables 

A = 02 - 0i and S = (f) 2 + <j>\ (7) 
and applies an expansion for small I 

A = Aq + ZAx + O^ 2 ), (8) 
£ = So + /S! + C(/ 2 ). (9) 

This analytical approach was complemented by a numerical investigation exploiting the 
Personal Supercomputer Circuit ANalyzer program PSCAN Jill, U We will not go into 



the details of calculation here, but concentrate on physical results. 

For vanishing external flux the elements behave in the same way as two weakly coupled 
elements; phase shift between voltages vanishes and elements oscillate oppositely with 
voltages 

vi = v 2 = — % — , (10) 

i + cos(Cos - d ) 

V k 



Co = \Ao - !> = const, v k = * ■ (11) 

The behavior is quite plausible, because in this case there is no current flow through the 
shunt. 

If the external flux does not vanish, the behavior is more complicated. To first order 
in I, time dependent voltage is obtained as follows: 

Co 2 _ , l Cosin(y/2) _ x 

i + cos((p/2) cos( s 2(i + cos(<^/2) cos Co 5 ) 2 

• / /0 x • 7 (■ . / /0 n , Co , io + cos(<^/2)cosCos\ 

sin <^/2 sin Cos «o + cos(^/2) + — — In — -— 

\ cos(y2/2) z + cos(v?/2) J 

=FCo(cos(^/2) + i cos Cos) 

where 



Co = sfi 2 - cos 2 (y./2). (13) 

This solution shows several specific features, (i) As in the case of weak coupling there 
exists a solution exhibiting phase locking for all values of the external field, (ii) For 
extremely strong coupling / — > both elements behave similar to free junctions, but with 
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the current i$ substituted by a kind of "effective current" depending on external flux, 
io/ cos((^/2). (ii) There is a phase shift between voltage phases being caused by the last 
term of Eq. fll2l) . However, contrary to the weak coupling case, this phase shift is very 
small over a wide range of external flux and only in the vicinity of the value ip = % the 
phase shift jumps to the value 5 = it. Within the range of validity of our perturbation 
theory it is possible to derive an explicit analytical expression for the phase shift 5: 



cos 5 



(b 2 + aia 2 ) 



b" + b 2 (al + al)+a\a 2 2 



(14) 



where oi, a 2 and b are the leading Fourier coefficients of the voltages fll2| 



Ol/2 



— Co _,_ , . / /0 x Co 
- 2 \ r— -j- T / sm(y/2) , 
\ «o + Co «o + Co 

Zsin%/2)(^ 

\ ?,nfn 



^oCo + Co \ 



io ~ Co 
*o + Co 



Co io -Cc 

4«g(i + Co) Vo + Cc 



(15) 
(16) 



Fig. 3 shows a plot of phase shift between junctions against external field. 

This can be compared with results from numerical simulation (including a small shunt 
capacitance (3c = 0.01; cf. Fig. 4). 

The limiting case (p — 7T is of special interest: 

vi = i - ysin(i s - S ) + ^sin2(z s - 5 ), (17) 
v 2 = i + y sm(i s -5 )+ l - sin2(i s - 5 ). (18) 
This clearly indicates the phase difference between both voltages to first order being equal 

to S = 7T. 

Comparison of Figs. 3 and 4 shows that even for I fh 1 where our analytical approx- 
imation is not longer valid results are not far from that of numerical simulation. Both 
figures show that already in this case the behavior is qualitatively different from what has 
been known for weakly coupled elements (dotted lines in Figs. 3 and 4). 



Within the interval < (f < tc Eq. fll4|) can be rewritten as 

5((p) = arctan(a 2 /&) — arctan(ai/6). (19) 
This shows, that for I < 1, but finite, the solution can be approximated by a step function 

S w n6(ip - (p*) (20) 

with the step at 

ip* = n-i l (21) 

(cf. Fig. 5). 
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Clearly ip* — > for I — > 0. This approximation might be useful considering more 
complicated arrays. 

To conclude, our investigation of two strongly coupled Josephson junctions has shown, 
that for every value of the external field there exists a regime showing phase-locking 
between voltage phases. Contrary to the case of weak coupling via large inductances, 
phase shift between voltages is negligibly small for a wide range of external flux for small 
inductive shunts, but shows a narrow peak centered around the value ip = n. If larger 
arrays will show a similar behavior, which has to be investigated in more detail, this would 
make them nearly insensible to parasitic external flux. 
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Figure Captions 



Fig 1. Circuit of two Josephson junctions with superconducting inductive SQUID coupling. 



Fig 2. Mean voltage phase shift 5 against normalized external field <p from analytical approxi- 
mation for weak inductive coupling (zq = 1-5). 



Fig 3. Mean phase shift 5 against normalized external filed <p for strong inductive coupling 
I = 0.1 and medium inductive coupling I = 1.0 obtained from analytical approximation (14) 
(to = 1.5). 



Fig 4. Mean phase shift 5 against normalized external filed cp for strong inductive coupling 
I = 0.1 and medium inductive coupling Z = 1.0 obtained by numerical simulation (iq = 1.5). 



Fig 5. Mean phase shift 5 against normalized external field <p for extremely strong coupling 
(0.001 < I < 0.1) obtained from analytical approximation ( |14"D (iq = 1.5). 
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Fig.l 



L 



Rn Rn 

[Si rSi 



/n 



in 



7 



Fig.2 
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Fig.3 
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Fig.4 
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Fig.5 
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